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I. INTRODUCTION 

Most theoretical and experimental studies of exclusive 
B — > X u £v decays have focused on B — > tt£v decays, 
while B — > plv and B + — > ui£ + v [l[ decays involving the 
vector mesons p and to have received less attention. Here 
£ is an electron or muon, and X refers to a hadronic 
state, with the subscript cor a signifying whether the 
state carries charm or is charmless. Measurements of 
the branching fraction of B — > p£v are impacted by an 
irreducible B — > X u lv background, typically the dom- 
inant source of systematic uncertainty. In studies of 
B + — > u)( + v that background can be suppressed to a 
larger degree, since the ui width is about 15 times smaller 
than that of the p. Extractions of the CKM matrix ele- 
ment \V u b\ from B + — > oj£ + v and B —J- p£v decay rates 
have greater uncertainties than those from B — > ixlv, 
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due to higher backgrounds and more complex form-factor 
dependencies. The persistent discrepancy between \y u b\ 
measurements based on inclusive and exclusive charmless 
decays is a motivation for the study of different exclusive 
B -> X u £v decays 0,01 . 

Measurements of B(B + — > lo£ + v) have been reported 
by Belle [1, H| ; a measurement by BABAR has been per- 
formed on a partial dataset jf|. In this analysis we use 
the full BABAR dataset to measure the total branching 
fraction B(B + — ¥ u£ + u) and partial branching fractions 
AB(B + —> lj£ + v) I Aq 2 in five q 2 intervals, where q 2 refers 
to the momentum transfer squared to the lepton system. 
The differential decay rate for B + —> loI + v is given by 



dr(s H 



dq 2 



= \v ub 



\H 4 



96n 3 m 2 B Cy 



(1) 



where p u is the magnitude of the to momentum in the 
B rest frame, mg is the B mass, and Gp is the Fermi 
coupling constant. The isospin factor cy is equal to \[2 
for B + lj£ + v §. As described in a related BABAR, 
paper the three helicity functions H^, H + , and H_ 
can be expressed in terms of two axial vector form factors 
A\ and and one vector form factor V, which describe 
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strong interaction effects, 



H±{q 2 ) = (m B + m u ) 



2m B p u 



H (q 2 ) = 



m B + m w 
Irn^xfq 1 



(m 2 B - ml - q^A^q 2 ) 
4m| pi 



(m B + m w ) : 



:A 2 (q' 



We compare the measured q 2 dependence of the decay 
rate with form factor predictions based on light-cone sum 
rules (LCSR) and the ISGW2 quark model We 
also use these form factor calculations and the measured 
branching fraction to extract | V u {, | . 



II. DETECTOR, DATA SET, AND SIMULATION 

The data used in this analysis were recorded with the 
BABAR detector at the PEP-II e + e~ collider operating at 
the T(45) resonance. We use a data sample of 426 fb -1 , 
corresponding to (467 ± 5) million produced BB pairs. 
In addition, we use 44 fb" 1 of data collected 40 MeV 
below the BB production threshold. This off-resonance 
sample is used to validate the simulation of the non- BB 
contributions whose principal source is e + e~ annihilation 
to qq pairs, where q = u,d, s, c. 

The PEP-II collider and BABAR detector have been de- 
scribed in detail elsewhere [ll|. Charged particles are 
reconstructed in a five-layer silicon tracker positioned 
close to the beam pipe and a forty-layer drift chamber. 
Particles of different masses are distinguished by their 
ionization energy loss in the tracking devices and by a 
ring-imaging Cerenkov detector. Electromagnetic show- 
ers from electrons and photons are measured in a finely 
segmented CsI(Tl) calorimeter. These detector compo- 
nents are embedded in a 1.5 T magnetic field of a su- 
perconducting solenoid; its steel flux return is segmented 
and instrumented with planar resistive plate chambers 
and limited streamer tubes to detect muons that pene- 
trate the magnet coil and steel. 

We use Monte Carlo (MC) techniques [HI, HH to simu- 
late the production and decay of BB and qq pairs and the 
detector response [lH , to estimate signal and background 
efficiencies and resolutions, and to extract the expected 
signal and background distributions. The size of the sim- 
ulated sample of generic BB events exceeds the BB data 
sample by about a factor of three, while the MC samples 
for inclusive and exclusive B — > X u lv decays exceed the 
data samples by factors of 15 or more. The MC sample 
for qq events is about twice the size of the qq contribution 
in the T(4S) data. 

The MC simulation of scmilcptonic decays uses the 
same models as in a recent BABAR analysis Q . The sim- 
ulation of inclusive charmless semileptonic decays B — > 
X u lv is based on predictions of a heavy quark expan- 
sion [l5| for the differential decay rates. For the simu- 
lation of B — > ■nlv decays we use the ansatz of [l6[ for 



the q 2 dependence, with the single parameter a B K set 
to the value determined in a previous BABAR analysis 
[l7| . All other exclusive charmless semileptonic decays 
B — > X u £u, including the signal, are generated with form 
factors determined by LCSR @,GJ]. For B Dlv and 
B — > D*lv decays we use parameterizations of the form 
factors [l|| [2(| based on heavy quark effective theory; 
for the generation of the decays B —t D**lv, we use the 
ISGW2 model [nj. 



III. CANDIDATE SELECTION 

In the following, we describe the selection and kine- 
matic reconstruction of signal candidates, the definition 
of the various background classes, and the application of 
neural networks to further suppress these backgrounds. 

The primary challenge in studying charmless semilep- 
tonic B decays is to separate signal decays from Cabibbo- 
favored B — > X c lv decays, which have a branching frac- 
tion approximately 50 times larger than that of B — > 
X u £v. A significant background also arises due to multi- 
hadron continuum events. 

Based on the origin of the candidate lepton we dis- 
tinguish three categories of events: 1) Signal candidates 
with a charged lepton from a true B + — > ojI + v de- 
cay; 2) BB background with a charged lepton from all 
non-signal BB events; 3) Continuum background from 
e + e~ — > qq events. The oj meson is reconstructed in its 
dominant decay, oj — > 7r + 7r~7r°. For each of the three cat- 
egories of events we distinguish correctly reconstructed 
oj — > 7T + 7r~7r° decays (true-w) from combinatorial-a; can- 
didates, for which at least one of the reconstructed pions 
originates from a particle other than the oj. 



A. Preselection 

Signal candidates are selected from events with at least 
four charged tracks, since a B + — > oj£ + v decay leaves 
three tracks and the second B in the event is expected to 
produce at least one track. The magnitude of the sum of 
the charges of all reconstructed tracks is required to be 
less than two, helping to reject events with at least two 
undetected particles. 

The preselection places requirements on the recon- 
structed lepton, oj meson, and neutrino from the B + — > 
ojI + v decay. At this stage in the analysis we allow for 
more than one candidate per event. 

The lepton is identified as cither an electron or muon. 
The electron identification efficiency is greater than 90% 
and constant as a function of momentum above 1 GeV, 
while the muon identification efficiency is between 65%- 
75% for momenta of 1.5-3 GeV. The pion misidentifi- 
cation rates are about 0.1% for the electron selector and 
1% for the muon selector. The lepton is required to have 
a momentum in the center-of-mass (cm.) frame greater 
than 1.6 GeV. This requirement significantly reduces 
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the background from hadrons that are misidcntified as 
lcptons, and also removes a large fraction of true lep- 
tons from secondary decays or photon conversions and 
from B — > X c iv decays. The acceptance of the detector 
for leptons covers momentum polar angles in the range 
0.41 < 9 < 2.54 rad. 

For the reconstruction of the decay oj — > Tr + ir~ir°, we 
require that the candidate charged pions are not iden- 
tified as leptons or kaons. The reconstructed oj mass 
must be in the range 680 < < 860 McV , and the 

7T° candidate is required to have an invariant mass of 
115 < m 77 < 150 McV. To reduce combinatorial w back- 
ground, we require minimum momenta for the three pion 
candidates, p^± > 200 MeV and p^a > 400 MeV, and 
also energies of at least 80 MeV for photons from the tt° 
candidate. 

The charged lepton candidate is combined with an oj 
candidate to form a so-called Y candidate. The charged 
tracks associated with the Y candidate are fitted to a 
common vertex Y vtx . This vertex fit must yield a y 2 prob- 
ability Prob(x 2 , Y vtx ) > 0.1. To further reduce back- 
grounds without significant signal losses, we impose two- 
dimensional restrictions on the momenta of the lepton 
and w. Each Y candidate must satisfy at least one 
of the following conditions on the cm. momentum of 
the lepton and oj: p* > 1.3 GeV, or p\ > 2.0 GeV, or 
Pi+plj > 2.65 GeV, where quantities with an asterisk 
refer to the cm. frame. These requirements reject back- 
ground candidates that are inconsistent with the phase 
space of the signal decay. The condition | cos#by| < 1.0, 
where cos6» B y = (2E* B E^ - M% - M^)/(2p* B p^) is the 
cosine of the angle between the momentum vectors of 
the B meson and the Y candidate, should be fulfilled 
for a well-reconstructed Y candidate originating from 
a signal decay [211 ] . The energy E* B and momentum 
p* B of the B meson are not measured event by event. 
Specifically, E* B ~ y/s/2, where yfs is the cm. en- 
ergy of the colliding beams, and the B momentum is 
derived as p* B = E* B — r?7^j . To allow for the finite 
resolution of the detector, we impose the requirement 
-1.2 < cos6>by < 1-1. 

The neutrino four-momentum is inferred from the 
missing energy and momentum of the whole event: 
(SmissjPmiss^ = (E e + e - , p e + e - ) - E u J2iPi), where 
E e + e - and p e + e - are the energy and momentum of the 
colliding beam particles, and the sums are performed over 
all tracks and all calorimeter clusters without an associ- 
ated track. If all tracks and clusters in an event are 
well-measured, and there are no undetected particles be- 
sides a single neutrino, then the measured distribution of 



the missing mass squared, m 



2 

miss 



E, 



Pmiss: P eaks 



at zero. We require the reconstructed neutrino mass 
to be consistent with zero, \m^ iss / (2E T[Liss )\ < 2.5 GeV, 
and the missing momentum to exceed 0.5 GeV. The 
polar angle of the missing momentum vector is also re- 
quired to pass through the fiducial region of the detector, 
0.3 < 6l mi ss < 2.2 rad. 

Other restrictions arc applied to suppress qq back- 



ground, which has a two-jet topology in contrast to BB 
events with a more uniform angular distribution of the 
tracks and clusters. Events must have i?2 < 0.5, where 
i?2 is the second normalized Fox- Wolfram moment [22j], 
determined from all charged and neutral particles in the 
event. We also require cos A# t hrust < 0.9, where A# t hrust 
is the angle between the thrust axis of the Y candidate's 
decay particles and the thrust axis of all other detected 
particles in the event. We require L2 < 3.0 GeV, with 
L2 = 'Yl li P*cos 2 0* , where the sum runs over all tracks 
in the event excluding the Y candidate, and p* and 9* 
refer to the cm. momenta and the angles measured with 
respect to the thrust axis of the Y candidate. 

We reject candidates that have a charged lepton and 
a low-momentum charged pion consistent with a B° — > 



D*~i+v, D*~ -> D { 'tt 



slow 



decay as described in [22 



The kinematic consistency of the candidate decay with 
a signal B decay is ascertained by restrictions on two 
variables, the beam-energy substituted B mass toes, 
and the difference between the reconstructed and ex- 
pected energy of the B candidate AE. In the lab- 
oratory frame these variables are defined as toes = 

y/( S /2 + p B ■ Pe+e- ) V^V - P| and AE = (P e+e - • 

P B - s/2)/,/s, where P B = (E b ,Pb) and P e + e - = 
(E e + e - ,p e +e- ) are the four- momenta of the B meson and 
the colliding beams, respectively For correctly recon- 
structed signal B decays, the AE distribution is centered 
at zero, and the toes distribution peaks at the B mass. 
We restrict candidates to -0.95 < AE < 0.95 GeV and 
5.095 < toes < 5.295 GeV. 



B. Neural Network Selection 

To separate signal candidates from the remaining back- 
ground we employ two separate neural networks (NN), 
to suppress qq background and B — > X c iv background. 
The qq NN is trained on a sample passing the preselec- 
tion criteria, while the B — > X c iv NN is trained on a 
sample passing both the preselection and the qq neural 
network criteria. The training is performed with signal 
and background MC samples. These NN are multilayer 
perceptrons that have two hidden layers with seven and 
three nodes. 

The variables used as inputs to the qq NN are i?2, L2, 
cosA6»thrust, cos6»by, m miss /(2£: miss ) , Prob(x 2 , V v tx), the 
polar angle of the missing momentum vector in the lab- 
oratory frame, and the Dalitz plot amplitude ^Daiitz = 
ot\p v + x p n -\, with the ir + and 7r~ momenta measured 
in the uj rest frame and scaled by a normalization fac- 
tor a. True u mesons typically have larger values of 
^Dalitz than combinatorial uj candidates reconstructed 
from unrelated pions. The B — > X c iv NN uses the same 
variables, except for cos A6> tm -ust, which is replaced by 
cos 6*n^, the hclicity angle of the lepton, defined as the 
angle between the momentum of the lepton in the rest 
frame of the virtual W and the momentum of the W in 
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the rest frame of the B. The data and MC simulation 
agree well for the NN input variables at each stage of the 
selection. The NN discriminators are chosen by max- 
imizing y^si„ + (1 — ?7bkg) 2 ! where e s ; g is the efficiency 

of the signal and f?bkg is the fraction of the background 
misidentihed as signal. 

The selection efficiencies for the various stages of the 
candidate selection for the signal and background com- 
ponents are given in Tabic |TJ After the preselection and 
NN selection, 21% of events in data contribute multi- 
ple B + — > loI + v candidates. The candidate with the 
largest value of Prob(x 2 , Y v t x ) is retained. For the re- 
maining candidates, the reconstructed 3-pion mass is re- 
quired to be consistent with the u> nominal mass [24j . 
|m 3Tr — m u \ < 23MeV. The overall signal efficiency is 
0.73% if the reconstructed candidate includes a true u> 
and 0.21% if it includes a combinatorial lo. The efficien- 
cies of the BB and qq backgrounds are suppressed by 
several orders of magnitude relative to the signal. 

TABLE I: Successive efficiencies (in %) predicted by MC 
simulation for each stage of the selection, for true- and 
combinatorial-ui signal, and backgrounds from BB and qq 
events. 



Source 


true-u; 
signal 


comb. -a; 
signal 


BB 


Qq 


Preselection 


1.9 


4.8 


0.0094 


0.00073 


Neural nets 


43 


17 


7.9 


11 


3-pion mass 


88 


26 


24 


30 


Total (product) 


0.73 


0.21 


0.00018 


0.000024 



C. Data-MC Comparisons 

The determination of the number of signal events relies 
heavily on the MC simulation to correctly describe the 
efficiencies and resolutions, as well as the distributions 
for signal and background sources. Therefore a signifi- 
cant effort has been devoted to detailed comparisons of 
data and MC distributions, for samples that have been 
selected to enhance a given source of background. 

Specifically, we have studied the MC simulation of the 
neutrino reconstruction for a control sample of B — > 
D*~i + v decays, with D*~ -> -D°tt s 7 ow and D° -> 
K + ir~ir°. This final state is similar to that of the 
B + —> lqI + v decay, except for the addition of the slow 
pion 7r~ and the substitution of a K + for a ir + . This con- 
trol sample constitutes a high-statistics and high-purity 
sample on which to test the neutrino reconstruction. We 
compare data and MC distributions for the control sam- 
ple and find good agreement for the variables used in 
the preselection and as inputs to the NN. We have also 
used this sample to study the resolution of the neutrino 
reconstruction and its impact on q 2 , toeSj an d AE. 




FIG. 1: Distribution of AE versus tties for true-u; signal 
MC. The 20 bins into which the plane is divided for the fit 
histogram are overlaid. 



IV. SIGNAL EXTRACTION 
A. Fit Method 

We determine the signal yields by performing an ex- 
tended binned maximum-likelihood fit to the observed 
three-dimensional AE-rriEs-q 2 distributions. The fit 
technique [25| accounts for the statistical fluctuations of 
the data and MC samples. 

For this fit the AE-tiies plane is divided into 20 bins, 
as shown in Fig. [TJ and the data are further subdivided 
into five bins in q 2 , chosen to contain roughly equal num- 
bers of signal events. The q 2 resolution is dominated 
by the neutrino reconstruction. It can be improved by 
substituting the missing energy with the magnitude of 
the missing momentum and by rescaling p m i ss to force 
AE = 0, q 2 m . r = [(Ei, ft) + S ■ (pmiss, Pmiss)] 2 , where 
(5 = 1 — AE / E m - lss . This correction to q 2 is used in the 
fit. 

We describe the measured AE-ttibs-Q 2 distribution as 
a sum of four contributions: B + — > ujt^v signal (both 
true-u; and combinatoric-u;), true-u; BB, true-u; qq, and 
the sum of the combinatorial-w background from BB and 
qq events. 

While the AE-ities shapes for the signal and true- 
u BB and qq sources are taken from MC samples, we 
choose to represent the dominant combinatorial-cj back- 
ground by the distributions of data events in the 
sidebands, thereby reducing the dependence on MC sim- 
ulation of these backgrounds. The normalization of these 
background data is taken from a fit to the 3-7T mass dis- 
tribution in the range 0.680 < m 37r < 0.880 GeV. To ob- 
tain a sample corresponding to the combinatorial-w back- 
ground from BB and qq events only, we subtract the MC 
simulated 7713^ contribution of the small combinatorial-u; 
B + — > Lot + v signal sample. To the resulting dis- 
tribution we fit the sum of a relativistic Breit-Wigncr 
convolved with a normalized Gaussian function, and the 
combinatorial background described by a second degree 
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polynomial. The resulting fit to the distribution for 
the all-g 2 sample is shown in Fig. [2] The \ 2 P er num- 
ber of degrees of freedom (dof) for the fits are within 
the range expected for good fits. The fitted background 
function is used to determine the weights to apply to the 
upper and lower sidebands to scale them to the expected 
yield of combinatorial-a; BB and qq background in the 
7713^ peak region. 
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FIG. 2: Fit to the distribution of mz-K for data from the 
all-g 2 sample, with MC combinatorial-o; signal subtracted. 
The dashed (red) and dotted (blue) curves describe the fit- 
ted peaking and combinatorial background functions, respec- 
tively, and the solid (black) curve is their sum. The peak and 
sideband regions are also indicated. 



The peak and two sideband regions are chosen to 
have a width of 46 MeV and are separated by 23 MeV, 
as indicated in Fig. [2j Since the normalization of the 
combinatorial-^ signal contribution depends on the fitted 
signal yield, which is a priori unknown, this component 
is determined iteratively. 

The fit has seven free parameters, five for the signal 
yields in each q 2 bin, and one each for the yields of the 
truc-w BB and qq backgrounds, the shapes of the dis- 
tributions are taken from MC simulations. The fitted 
yields are expressed as scale factors relative to the de- 
fault yields of the MC simulation. The total signal yield 
is taken as the sum of the fitted yields in the individual 
q 2 bins, taking into account correlations. 

B. Fit Results 

The fitting procedure has been validated on pseudo- 
experiments generated from the MC distributions. Wc 
find no biases and the uncertainties follow the expected 
statistical distribution. 

The yields of the signal, true-cj BB, and true-cj qq com- 
ponents obtained from the binned maximum-likelihood 
fit to AE-m-Es-q 2 are presented in Table [TTJ Projections 
of the fitted distributions of tties f° r the ah-<7 2 fit and for 
the five q 2 bins fit are shown in Fig. [31 The agreement 



between the data and fitted MC samples is reasonable 
for distributions of AE, tties, an d q 2 , as indicated by the 
X 2 /dof of the fit, 106/93, which has a probability of 16%. 
The fixed combinatorial-^ background yield accounts for 
83% of all backgrounds. The correlations among the pa- 
rameters are listed in Table UTil The strongest correlation 
is —72%, between the signal and qq yields in the first q 2 
bin, which contains most of the qq background. The cor- 
relation between signal and BB background is strongest 
in the last q 2 bin, —40%, because of a large contribution 
from other B — > X u lv decays. Correlations among signal 
yields are significantly smaller. 

The branching fraction, B(B + — > oj£ + is), averaged 
over electron and muon channels, is defined as B(B + —> 
ut+v) = Ei(Nf s /e( B )/(4f±N BE ), where A/f g refers to 
the number of reconstructed electron and muon signal 
events in q 2 bin i, ef 6 is the reconstruction efficiency, f± 
is the fraction of B + B~ decays in all BB events, and 
iVgg is the number of produced BB events. The factor 
of 4 comes from the fact that B is quoted as the average 
of i = e and [i samples, not the sum, and the fact that 
cither of the two B mesons in the B + B~ event may decay 
into the signal mode. The q 2 resolution in the signal re- 
gion is 0.36 GeV 2 , smaller than the width of the q 2 bins. 
To account for the finite q 2 resolution, the background- 
subtracted, efficiency-corrected spectrum is adjusted by 
deriving from the signal MC the ratio of the true and 
reconstructed q 2 spectra, (dB / dg 2 ruo ) / (dB / dq 2 cco ) . This 
ratio differs by 9% at low q 2 and considerably less at 
higher q 2 . The partial and total branching fractions listed 
in Table IIVI are corrected for the effects of finite q 2 reso- 
lution and efficiency. 



V. SYSTEMATIC UNCERTAINTIES 

Table [V] summarizes the contributions to the system- 
atic uncertainty. The event reconstruction systematic 
uncertainties are most sensitive to the neutrino recon- 
struction, which depends on the detection of all of the 
particles in the event. To assess the impact of the un- 
certainty of the measured efficiencies for charged tracks, 
the MC signal and background samples are reprocessed 
and the analysis is repeated, after tracks have been elim- 
inated at random with a probability determined by the 
uncertainty in the tracking efficiency. Similarly, we eval- 
uate the impact from uncertainties in the photon recon- 
struction efficiency by eliminating photons at random as 
a function of the photon energy. Since a K® leaves no 
track and deposits only a small fraction of its energy 
in the calorimeter, the reconstruction of the neutrino is 
impacted. The uncertainty on the K® MC simulation 
involves the shower energy deposited by the K® in the 
calorimeter, the K® detection efficiency, and the inclusive 
K® production rate as a function of momentum from BB 
events. 

The impact of the changes to the simulated background 
distributions which enter the fit arc smaller than for 
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TABLE II: Number of events and their statistical uncertainties, as determined from the fit, compared with the number of 
observed events in data. The combinatorial-^; background yields are fixed in the fit; the quoted uncertainties are derived from 
the sideband subtraction. 

q 2 range (GeV 2 ) 0-4 4-8 8-10 10-12 12-21 0-21 

All signal 257 ± 72 238 ± 44 161 ± 32 177 ± 32 293 ± 57 1125 ± 131 

True-w signal 238 209 136 137 168 869 

Comb.-w signal 19 28 25 40 125 256 

BB (true-w) 105 ± 19 192 ± 34 154 ± 27 195 ± 34 411 ± 73 1057 ± 187 

qq (true-w) 409 ± 96 145 ± 34 65 ± 15 34 ± 8 64 ± 15 716 ± 167 

Comb.-w bkgd. 1741 ± 23 1818 ± 24 1240 ± 20 1520 ± 22 3913 ± 35 10232 ± 57 

Data 2504 ± 50 2433 ± 49 1605 ± 40 1858 ± 43 4738 ± 69 13138 ± 115 



< q 2 < 4GeV 2 
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FIG. 3: Distributions of tties after the fit and the ratio of the data to the fitted predictions, for five separate q 2 bins and the 
full q 2 range, in the AE signal band, —0.25 < AE < 0.25 GeV. The points represent data with statistical uncertainties, while 
the stacked histograms represent the sum of fitted source components, signal (white), true-o; BB (light gray), true-o; qq (dark 
gray), and combinatorial-tj background (diagonally thatched). 



the signal, since the large combinatorial backgrounds are 
taken from data, rather than MC simulations. As an es- 
timate of the impact of these variations of the MC sim- 
ulated distributions on the q 2 dependent signal yield, we 
combine the observed reduction in the signal distribu- 
tion with the impact of the changes to qq and BB back- 
grounds on the signal yield, taking into account the cor- 
relations obtained from the fit (see Table III). Since the 
correlations between signal and backgrounds are small at 



high q 2 , the impact of the uncertainties in the background 
are also modest. This procedure avoids large statistical 
fluctuations of the fit procedure that have been observed 
to be larger than the changes in the detection efficien- 
cies. However, this procedure does not account for the 
small changes in the shape of the distributions, and we 
therefore sum the magnitude of the changes for signal 
and background, rather than adding them in quadrature 
or taking into account the signs of the correlations of the 
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TABLE III: Correlations among the fit scale factors pi for 
the simulated source s and q 2 bin k. The scale factors for qq 
and BB apply to the full q 2 range. 



P' 



P 



P _ 
P BS 
pf" 
pf" 

P3 

pt" 

P5 



1.000 



-0.466 
1.000 



uitv 

Pi 


P2 


uitv 
P3 


Pi 


ujiv 
P5 


0.724 


-0.106 


-0.031 


0.051 


0.088 


0.223 


-0.249 


-0.253 


-0.284 


-0.401 


1.000 


0.121 


0.061 


0.001 


-0.011 




1.000 


0.105 


0.094 


0.128 






1.000 


0.088 
1.000 


0.121 
0.125 
1.000 



TABLE IV: Measured B + — >• uil + v branching fraction and 
partial branching fractions in bins of q 2 with statistical and 
systematic uncertainties. 



(f (GeV) AB (x!0~ 4 ) 

0-4 0.214 ± 0.060 ± 0.024 

4-8 0.200 ± 0.037 ± 0.010 

8-10 0.147 ± 0.029 ± 0.010 

10-12 0.169 ± 0.031 ± 0.010 

12-21 0.482 ± 0.093 ± 0.038 

0-12 0.730 ± 0.083 ± 0.054 

0-21 1.212 ± 0.140 ± 0.084 



the combinatorial-w signal is not subtracted from the 
TO 37r sidebands, and the signal fit parameter is set to 
scale only the true-o; signal contribution. Twenty per- 
cent of the difference between the nominal and uncor- 
rected results is taken as the systematic uncertainty; it 
is largest for 12 < q 2 < 21 GcV 2 because the fraction 
of combinatorial-w signal in this q 2 bin is large. The 
sideband event yields determined from the fit are 
varied within their fit errors to determine the statisti- 
cal uncertainty on the combinatorial-^ background. The 
uncertainty in the chosen 7713^ ansatz is assessed by re- 
peating the fits, replacing the nominal functions for 
the peak and background components. For the back- 
ground component, we use a third instead of a second 
degree polynomial. For the peaking component, we use a 
Gaussian function in place of a relativistic Breit-Wigncr 
convoluted with a Gaussian function. The systematic er- 
ror from the ansatz is taken as the sum in quadrature 
of the change in signal yield for each of these functional 
variations. 



The branching fraction depends inversely on the value 
of N B g, which is determined with a precision of 1.1% [27j |. 
At the T(4S) resonance, the fraction of B + B~ events is 
measured to be f± = 0.516 ± 0.006 [24| . with an uncer- 
tainty of 1.2%. 



signal and backgrounds in a given q 2 bin. 

We assign an uncertainty on the identification effi- 
ciency of electrons and muons, as well as on the lepton 
and kaon vetoes of the u> daughter pions, based on the 
change in signal yield after varying the selector efficien- 
cies within their uncertainties. 

The uncertainty in the calculation of the LCSR form 
factors impacts the uncertainty on the branching fraction 
because it affects the predicted q 2 distribution of the sig- 
nal and thereby the fitted signal yield. We assess the im- 
pact by varying the form factors within their uncertain- 
ties. We include the uncertainty on the branching frac- 
tion of the u decay, B(u -> tt+tt-tt ) = (89.2±0.7) x 10~ 2 
24 1 . To evaluate the uncertainty from radiative correc- 
tions, candidates are reweightcd by 20% of the difference 
between the spectra with and without PHOTOS [26jj, 
which models the final state radiation of the decay. 

The uncertainty on the true-w backgrounds has a small 
impact on the signal yield since these components repre- 
sent a small fraction of the total sample. To assess the un- 
certainty of the AE-rriBS-Q 2 shapes of the true-w qq and 
true-w BB samples, the fit is repeated after the events 
are reweightcd to reproduce the inclusive uj momentum 
distribution measured in BB and qq events. We also as- 
sess the uncertainty on the modeling of the semileptonic 
backgrounds by varying the branching fractions and form 
factors of the exclusive and inclusive B — > X u lv 24j and 
B — > X r iv backgrounds Q within their uncertainties. 

To assess the uncertainties that result from the MC 
prediction of the m 37r distribution of the combinatorial- 
ui signal, we use the uncorrected distribution, in which 



TABLE V: Systematic uncertainties in % on the branching 
fraction. 



q 2 range (GeV 2 ) 


0-4 


4-8 


8-10 


10-12 


12-21 


0-21 


Event reconstruction 














Tracking efficiency 


3.9 


1.5 


2.8 


2.3 


1.1 


2.0 


Photon efficiency 


2.0 


1.7 


3.3 


1.1 


0.6 


1.5 


Kl prod. /interactions 


4.8 


1.8 


2.5 


1.1 


1.4 


1.9 


Lepton identification 


1.6 


1.5 


1.5 


1.2 


1.2 


1.3 


K/l veto of ui daughters 


1.7 


1.7 


1.7 


1.7 


1.8 


1.7 


Signal simulation 














Signal form factors 


6.3 


1.5 


1.1 


2.9 


4.6 


4.8 


B(tJ -»■ TT + TT-TT ) 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 


Radiative corrections 


0.4 


0.3 


0.2 


0.1 


0.2 


0.2 


True-o; background 














qq AE-niES-q 2 shapes 


2.6 


0.1 


0.4 


0.2 


0.3 


0.5 


BB AE-rriES-q 2 shapes 


2.0 


0.9 


1.8 


0.2 


0.1 


0.8 


B ->• X c lv B and FF 


0.2 


0.6 


0.3 


0.2 


0.2 


0.2 


B X u lv B and FF 


0.3 


0.4 


0.4 


0.3 


0.5 


0.4 


Comb.-u sources 














Signal 7713^ distribution 


0.6 


0.5 


0.4 


1.1 


3.7 


1.5 


Bkgd. yield, stat. error 


4.2 


1.0 


0.9 


0.9 


2.0 


1.7 


Bkgd. yield, ansatz error 


1.7 


2.2 


2.7 


2.7 


3.5 


0.9 


B production 














BB counting 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


f± 


1.2 


1.2 


1.2 


1.2 


1.2 


1.2 


Syst. uncertainty 


11.1 


5.2 


6.8 


5.8 


7.9 


6.9 


Stat, uncertainty 


28.1 


18.7 


20.0 


18.1 


19.4 


11.6 


Total uncertainty 


30.2 


19.4 


21.1 


19.0 


20.9 


13.5 
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VI. RESULTS AND CONCLUSIONS 

We have measured the branching fraction, 

B(B+ -> uj( + v) = (1.21 ± 0.14 ± 0.08) x 10~ 4 , (2) 

where the first error is statistical and the second is sys- 
tematic, based onll25± 131 observed signal candidates. 
Here, £ indicates the electron or muon decay mode and 
not the sum over them. The measured partial branching 
fractions are presented in Table ITvl and are compared to 
the predictions from two form factor calculations in Fig. 
21 These QCD predictions have been normalized to the 
measured branching fraction. 



,y 1 = / Ag (gmin>gmax) 
V T + A C(<?min 

AC(A„,0 = ^/f"^fFW. (3) 

Table Ell lists the values of AC and \V ub \ for LCSR 
and ISGW2 in different ranges of q 2 . LCSR calculations 
are more accurate at low q 2 , while ISGW2 predictions 
are more reliable at high q 2 . Both form factor calcula- 
tions arrive at very similar values for \V u b\. These values 
of | Kit, | are consistent with the more precisely measured 
values from B — > ■niv decays [29| . 



> 
d> 
O 

CM 

< 

< 




20 

q 2 (GeV 2 ) 



FIG. 4: Partial branching fractions (points with error bars) 
with respect to q 2 . The data are compared with the predic- 
tions from light-cone sum rules (LCSR) [i| and a quark-model 
calculation (ISGW2) [Tcj. The uncertainty band (shaded) is 
given for the LCSR calculation. 



Neglecting the theoretical uncertainties, the % 2 /NDF 
of the measured distribution relative to the LCSR pre- 
diction Q is 2.4/4, corresponding to a y 2 probability of 
67%; relative to the ISGW2 prediction jnj the X 2 /NDF 
is 4.2/4, with a \ 2 probability of 40%. Within the large 
experimental uncertainties, both the LCSR and ISGW2 
form factor calculations are consistent with the data. The 
uncertainties of the ISGW2 form factor calculation are 
not available. The uncertainties of the LCSR calcula- 
tion were estimated by the authors to vary linearly as 
a function of q 2 ; i.e., a dB / dq 2 /(dB/dq 2 ) = 21% + 3% x 
q 2 /(UGeV 2 ), for the B -> p£v decays It is assumed 
that this estimate is also valid for B + — > u£ + v decays. 

The value of \V u b\ can be determined from the mea- 
sured partial branching fraction, the B + lifetime r + = 
(1.638 ± 0.011) ps [HI, and the integral AC of the pre- 
dicted differential decay rate: 



TABLE VI: \V u b\, determined from two form factor calcula- 
tions of AC, in different ranges of q 2 . The first uncertainty 
is experimental (the sum in quadrature of statistical and sys- 
tematic); the second uncertainty is from theory, and is only 
available for LCSR. 



;GeV 2 ) AC(ps- 1 ) |K, 



fxHT 



LCSR [8] 



ISGW2 [10] 



0-12 
12-21 
0-21 

0-12 
12-21 
0-21 



3.9 ± 0.9 
3.2 ± 0.8 



3.37 ± 0.23 ± 0.38 
3.04 ± 0.32 ± 0.37 



7.1 ± 1.7 3.23 ± 0.22 ± 0.38 



3.6 
3.4 
7.0 



3.51 ± 0.24 
2.94 ± 0.31 
3.24 ± 0.22 



The value of B(B + — > ujI + v) measured in this analy- 
sis supersedes the previous BAB4R measurement @ based 
on a smaller data sample, and is in excellent agreement 
with a recent result [301 based on the full BABAR data 
set. The principal difference between this analysis and 
the previous ones is that the combinatorial-w background 
is taken from the sideband of the data distribution 
rather than from MC simulation. Although the domi- 
nant systematic uncertainties from event reconstruction 
cannot be avoided, this procedure substantially reduces 
the reliance on the MC simulation of this largest source 
of background. 

Currently, the QCD predictions of the form factors, 
and in particular their uncertainties, have limited preci- 
sion for B + — > uj£ + v and B — > piv decays. These form 
factor uncertainties impact |Kb| derived from B(B + —> 
uj£ + v). In the future, form factor calculations with re- 
duced uncertainties combined with improved branching 
fraction measurements would enable tests and discrim- 
ination among different predictions as a function of q 2 , 
and thereby improve the determination of |K&|- 
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